Correlated few-particle states in artificial bipolar molecule 
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We investigate the ground and excited states of a bipolar artificial molecule composed of two 
vertically coupled quantum dots containing different type of carriers - electrons and holes - in 
equilibrium. The approach based on exact diagonalization is used and reveals an intricate pattern 
of ground-state angular momentum switching and a rearrangement of approximate single-particle 
levels as a function of the inter-dot coupling strength. 
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In recent decades, artificially created zero-dimensional 
systems — the so-called quantum dots or artificial atoms 
- have been a major focus of active research in con- 
densed matter physics M. Most of the interest arises 
from our ability to control the addition or extraction of 
even a single particle fi ||] as well as from the fact that 
these artificial structures often offer a rich set of novel 
features unrivaled by those of their natural counterparts. 
For example, while the nature of atomic spectra is mostly 
governed by the quantization of motion of a single parti- 
cle in the field of the nucleus, in quantum dots the crucial 
role of the quantizing magnetic field and strong Coulomb 
interaction between the electrons leads to the formation 
of intriguing 'magic' states & ^ related to the compact 
states of non-interacting composite fermions Q || . 

The ability to create and study double-layer systems 
has added a new dimension to the physics of quantum 
dots. Two vertically stacked quantum dots make an ex- 
ample of an artificial molecule. The observation of a novel 
incompressible state induced by the inter-layer coupling 
H motivated a study of correlated few-particle states in 
vertical double-dot systems M. In this work, the inter- 
dot interactions were analysed in terms of an isospin- 
1/2, and a new class of incompressible ground states was 
found. 

In the present paper, we examine the correlated few- 
particle states of an artificial bipolar molecule consist- 
ing of two vertically stacked quantum dots populated by 
carriers of different kinds - electrons and holes, respec- 
tively. Such a system could be structured, for example, in 
bilayer-bipolar heterostructures containing electron and 
hole layers in equilibrium. These structures have been 
realized in the crossed gap InAs/GaSb system |q| as well 
as in biased GaAs/AlGaAs heterostructure ||, [HJ where 
electron and hole layers form on the opposite sides of an 
AlGaAs barrier. Most of the interest in such systems 
stems from the possibility (at least, in principle) of the 
formation of Bose-Einstein condensate of indirect exci- 
tons 0. 

While the superfluid state has not been demonstrated 
so far, a number of other interesting effects of electron- 
hole coupling and hybridization have been predicted and 



observed. For example, the Coulomb-drag experiments 
|[L0[ revealed the inter-layer momentum transfer rates 
up to an order of magnitude larger than in similar 
electron-electron systems, and magnetotransport calcula- 
tions |l3| established the presence of coupling-dependent 
current that remains large even when the Fermi energy 
lies in the hybridization gap. 

In the present paper, we consider two vertically cou- 
pled infinitesimally thin circularly-symmetric quantum 
dots placed into a strong perpendicular magnetic field. 
The confining potentials of the dots are assumed to be of 
a simple parabolic shape, and we choose the base fre- 
quencies of the two dots to be inversely proportional 
to the effective masses of the respective carriers, i.e. 
co e m e = UhTTih- This choice implies the equality of 
the quantum dot oscillator lengths Iq — ^Jh/m a uj a for 
a = e,h (a plausible assumption) as well as the equal- 
ity of renormalized oscillator lengths in magnetic fields 



of arbitrary strength since I 4 
cyclotron radius / 



+ and the 



\JhcjeB does not depend on the 
material parameters. The above constraint also consid- 
erably simplifies the calculation of Coulomb matrix ele- 
ments between the electron and hole states (see below). 

In the limit of strong magnetic fields, when the cy- 
clotron frequencies of the two types of carriers u> c / e \ and 
uj c (h) exceed the respective base dot frequencies (in prac- 
tice, B > 2-3 T) we can restrict the single-particle bases 
in the dots to the essentially one-dimensional sets of 
the lowest-Landau-level states. The energies of these 
states are spaced equidistantly by e a = 7io;„/<j c ( a ) with 
a = e,h. We note that due to the inverse proportionality 
of the cyclotron frequency to the effective mass we have 
e e /e h = m h /m e . 

These states are labelled by a single index n = 0, 1, 2 . . . 
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where the upper sign is to be used for electrons and the 
lower one for holes. That is, the electrons can only have 
negative or zero angular momenta and vice versa. More- 
over, since the particles in strong magnetic fields are spin- 
polarized we can neglect the presence of the spin degree 



2 



of freedom and also drop the constant Zeeman energy. 

For a many particle system, a convenient basis is given 
by the set of all possible Slater determinants (equiva- 
lently, ordered strings of the associated creation opera- 
tors) composed of the functions (|l]). Thus we will denote 
possible configurations of the molecule by enumerating 
the single-particle quantum numbers rij in ascending or- 
der, first for electrons and then, separated by a semicolon, 
for holes. E.g. the ground state of three non-interacting 
electrons and holes is given by {012;012}. Due to the 
cylindrical symmetry of the molecule the many-body 
states of different total angular momenta M t will not 
be mixed, thus we will be able to perform diagonaliza- 
tions separately in subspaces of a given M t . As a matter 
of fact, in our basis the single-body part of the Hamilto- 
nian is diagonal and the remaining task is the calculation 
and diagonalization of the Coulomb interactions. 

Before proceeding any further, it is convenient to intro- 
duce handy dimensionless units thus also clarifying the 
relevant energy and length scales of the problem. From 
now on we will measure all the energies in e e and the 
lengths in /. Then the spacing of the electronic single- 
particle is set equal to 1, and the spacing of the hole levels 
is m e /m/j. The interaction part of the Hamiltonian be- 
comes 
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Here, A represents the vertical separation between the 
two dots, and 
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expresses the relative strength of inter-particle interac- 
tions in terms of a ratio of the characteristic lengths, a* B 
being the effective electronic Bohr radius. 

We calculate the matrix elements of the Coulomb in- 
teraction between various many-body states by decom- 
posing them into a sum of two-particle Coulomb integrals 



J= Jefn y^r 2 C>i)C>2Wn 3 (ri)<K(r 2 ) (4) 

with v = |i"i — r2 1 — 1 if both particles belong to the same 
layer, and v — |ri — r 2 + A| _1 otherwise. In the former 
case the integral is easily evaluated analytically using the 
complex- variable technique due to Girvin and Jach Q to 
separate the CM and relative coordinates (see formula (7) 
in their paper). When the particles ri and r 2 belong to 
different layers the only modification is the replacement 
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in the integral over the radial part 
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FIG. 1: Angular-momentum resolved spectrum of an artificial 
molecule in the weak-coupling limit. The stable low-energy 
'magic' states are marked by open dots and joined by a dashed 
line. 



of the relative coordinate p which has to be computed 
numerically. 

We numerically calculated energy levels and many- 
particle states of neutral systems consisting of equal num- 
bers (1, 2 and 3) of electrons and holes in the dots. How- 
ever, in order to keep the presentation concise, we will 
concentrate on the most illustrative results obtained for a 
six-particle artificial molecule. Systems involving smaller 
number of particles show qualitatively similar behaviour. 

Weakly correlated dots. When the distance between 
the two dots is relatively large — such as A = 3.3 (in the 
units I) to be used in the present example — the inter- 
dot interactions between the particles are rather weak 
compared to the strong Coulomb repulsion within the 
individual dots which is decisive in the determination of 
the overall behaviour of the molecule. In Fig. [l] we plot 
the angular-momentum resolved spectrum of the system 
showing 15 lowest states for each value of M t at A = 4.0. 
One sees that for the values of M t equal to multiples of 
3 the lowest-energy states (marked by the open circles) 
are typically considerably lower than those correspond- 
ing to the neighbouring values of the angular momen- 
tum. This is a consequence of the fact that the individ- 
ual dots 'prefer' to be in the well-known 'magic' states. 
For three particles (electrons or holes) the predominant 
single-particle configurations in these stable many-body 
states are {012}, {123}, {234}, etc. It will be instructive 
to compare these configurations with their counterparts 
observed in strongly coupled dots below. 

Since the effective mass of the holes is considerably 
larger than that of electrons (we use mt — 6.7m e ) the 
Coulomb repulsion promotes them to states of higher an- 
gular momentum more easily. In Fig. |l| we see that the 
ground state has M t = 6, and the dashed lowest-energy 
curve rises much more sharply in the direction of negative 
angular momenta. 
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FIG. 2: The dependence of relative ground-state energies 
with Mt — 0, 3, 6, 9 on the interaction strength A. The inset 
shows the switching of the angular momentum of the absolute 
ground-state. 



FIG. 3: The evolution of the relative ground-states with —3 < 
Mt < 7 versus the coupling strength A. The ground state 
energy is subtracted from all the curves. 



While the angular momentum of the ground state is 
generally expected to grow with increasing Coulomb re- 
pulsion, this growth is not monotonous. The inset of 
Fig. |^ shows the evolution of the ground-state M t ver- 
sus A. Notice that the angular momentum reaches the 
value M t = 9 in three monotonous steps AM t — 3, and 
then switches back and forth to the value M t = 6. To 
understand this kind of behaviour we also investigated 
the evolution of the relative ground-states (i.e. lowest- 
energy states of a given M t ) of the involved angular mo- 
menta in Fig. 0. We see that the A-dependences of these 
energies typically consist of segments of nearly straight 
lines connected by abrupt breaks associated with the 
change of the nature of the states. For example, we 
observe two breaks in the M t — curve at A « 4 and 
A w 9. At A w 4 the molecule minimizes its energy 
by switching from the dominant configuration {012;012} 
into {123;123} followed by the changeover into {234;234} 
at A sa 9. The curves corresponding to higher values of 
Mt also show similar behaviour which in this particular 
case results in a close competition between the Mt = 6 
and Mt — 9 lowest-energy states and switching of the 
absolute ground-state in a rather extended range of A 
values. 

Strongly correlated dots. Bringing the two dots closer 
(we will consider the case of A = 0.2) results in an en- 
tirely different character of their energy spectra. The 
Coulomb attraction between particles of different kind 
wins over the intra-dot repulsion and leads to the forma- 
tion of rather strongly bound and weakly interacting ex- 
citons. The energies of the relative ground-states of given 
angular momenta decrease rapidly and monotonously 
with increasing A and no longer contain abrupt breaks. 
The absolute ground state always has the total angular 
momentum M t = 0. For the clarity of presentation, in 
Fig. pi we plot the differences of the lowest-energy states 



of M t = — 3 . . . 7 and the absolute ground state (M t = 0). 
In this figure we see a clear rearrangement of the struc- 
ture of the energy levels from the weakly interacting limit 
to the strongly-interacting limit (high A). At A — > the 
levels corresponding to the positive momenta are spaced 
more closely that the negative-momentum levels due to 
the higher hole mass. For sufficiently high A values the 
energy levels regroup so that the levels corresponding to 
the angular momenta Mt and — Mt form closely spaced 
pairs. 

A very interesting feature of the correlated electron 
and hole states in this regime can be observed by looking 
more closely at the dominant single particle configura- 
tions of the many-body states. In the table [| we compare 
the configurations corresponding to the relative ground 
states for M t — ... 7 in the weakly and strongly cou- 
pled dots. Since we are looking at the positive momen- 
tum states the electronic subsystem is in the energeti- 
cally most favourable state {012} of the total momentum 
M e = — 3 while the holes contribute the angular momen- 
tum Mh = M t + 3. However, we see that the pattern of 
filling of the hole states is very distinct in the two regimes. 
In the weak coupling limit, the holes favour the 'tradi- 
tional magic' filling scheme in which they occupy a com- 
pact set of adjacent angular momentum orbitals, such as 
{123} in the case of Mh = 6. In contrast, when the holes 
are in the Mh = 6 states in the strong coupling limit, the 
configuration {123} becomes the least favoured one. In- 
stead, the dominant configuration is {015} in which two 
of the holes occupy the orbitals of the lowest possible 
angular momentum while the third one is placed into a 
remote orbital n = 5. 

This new type of systematics of the single-particle con- 
figurations can be understood by realising that singling 
out one particle into a remote high-n orbital becomes 
energetically favourable in the case when the spacing be- 
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TABLE I: The predominant single-particle configurations of 
the relative ground states. 

M t Weak inter-dot coupling Strong inter-dot coupling 






{012;012} 


{012;012} 


1 


{012;013} 


{012;013} 


2 


{012;023} 


{012;014} 


3 


{012;123} 


{012;015} 


4 


{012;124} 


{012;016} 


5 


{012;134} 


{012;017} 


6 


{012;234} 


{012;018} 


7 


{012;235} 


{012;019} 





0123456 0123456 0123456 

FIG. 4: A schematic plot of single-particle energies versus an- 
gular momentum. Panel (a) displays the case of equidistant 
level spacing, (b) corresponds to the case of energy levels con- 
densing downwards, and (c) represents the case of upwards- 
condensing energy levels. Closed dots in panels (b) and (c) 
depict the lowest-energy configurations of total momentum 
M h = 6. 



tween energy levels is decreasing with increasing n, while 
in the opposite case the preferred configuration is the 
compact one. We illustrate this point in Fig. |J. Panel 
(a) shows the dependence of the single-particle energy on 
its angular momentum for non-interacting particles in a 
quantum dot. Since the energy levels are equidistant, all 
three ways of accommodating three particles of the total 
momentum M h = 6 — {015}, {024} and {123} — have 
equal total energies. However, if the energy levels are 
condensing downwards [panel (b)], the preferred config- 
uration is the compact {123} one as we observe in the 
weak-coupling limit. On the other hand, if the energy 
levels are condensing upwards, as shown in panel (c), the 



favoured configuration is {015} thus reproducing our re- 
sult in the strong-coupling limit. 

In conclusion, we have investigated the many-body 
states of a bipolar artificial molecule consisting of two 
vertically aligned quantum dots in two distinct regimes. 
In the weak-coupling limit we found an intricate pattern 
of switching of the ground state angular momentum as a 
function of a perpendicular magnetic field, while in the 
opposite strong-coupling limit we observe the formation 
of a new class of energetically favourable states. 
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